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FOREWORD 


This research was performed under Program Element 61102H, Project 5710, 
Subtask PA051, and was funded by the Defense Atomic Support Agency (DASA). 


Inclusive dates of research were July 1968 to December 1968. The report 
was submitted 6 October 1969 by the Air Force Weapons Laboratory Project Officer, 
Captain Robert E. Wiley (WLTH). 


The latitude-longitude data defining the continental outlines was obtained 
from the National Center for Atmospheric Research along with a sample program 
for plotting the modified mercator map. Subroutine ELIM is the only routine 
which remains essentially unchanged from the NCAR version. 


Information in this report is embargoed under the U.S. Export Control Act 
of 1949, administered by the Department of Commerce. This report may be released 
by departments or agencies of the U.S. Government to departments or agencies of 
foreign governments with which the United States has defense treaty commitments, 
subject to approval of AFWL (WLTH). 


This report has been reviewed and is approved. 
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ABSTRACT 


(Distribution Limitation Statement No. 2) 


WORLD is a general-use computer subroutine which draws maps containing 
continental outlines and overlays data on the map. Projection of maps 
may be modified mercator, equidistant polar, or hemispherical. 
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SECTION I 


INTRODUCTION 


Frequently one generates or obtains data which readily lends itself to plots 
that emphasize the positional relationship between the data and the earth. 
Generally this means that one must plot the data on a standard map of the earth 
(e.g., Calcomp paper W2), and be limited by the fixed scale, or draw a simple 
grid labeling the latitudes and longitudes and neglect the geographical features 
of continents, etc. Use of subroutine WORLD allows one to draw a map of reasona- 
ble size of any section of the earch (excluding the polar regions) and plot data 
on this map. The projection of the plot may be either modified mercator, equi- 
distant polar, or hemispherical. Hemisphere projections may be done as if viewed 


from any point above the earth. 
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SECTION II 


HOW TO USE THE ROUTINE 


The calling statement for WORLD is as follows: 
CALL WORLD (XLONG, YLAT, N, XSIZE, YSIZE, XEAST, XWEST, YNORTH, YSOUTH, MAP, 
IPLOT, IEND, IXAX, IX, IYAX, IY, ITLE, IT, LABEL, IL) where 


XLONG 


XEAST 


XWEST 


YNORTH 


YSOUTh 


The array containing the logitude data. The data must be in degrees, 
but may be in the range 0° to 360° or -180° to +180° (- indicates 
west longitude and + east) 


The array containing the latitude data. The data must be in degrees 
and in the range -90° to +90° (- indicates south latitudes and + 
north) 


The number of data points in the arrays XLONG and YLAT. If N= 0 
no data will be plotted, although any operations specified by other 
parameters will be performed. 


The length of the longitude axis in inches. 


The length of the latitude axis in inches. For polar or hemispheri- 
cal projection YSIZE is the diameter of the map. 


The longitude of the eastern boundary of the map. Units must be 
consistent with definition of XLONG. Not used for polar map. For 
MAP=5, longitude of viewing position. 


The longitude of the western boundary of the map. Units must be 
consistent with definition of XLONG. Not used for polar map. For 
MAP=5, latitude of viewing position. 


The northern boundary of the map. Units must be consistent with the 
definition of YLAT. Boundary of north polar map. 


The southern boundary of the map. Units must be consistent with the 
definition of YLAT. Boundary of south polar map. 
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IPLOT 


TEND 


IXAX 


IX 


IYAX 


IY 


ITLE 


Indicates the type of projection. The values which MAP may have 


and what each indicates are 

1. Modified mercator projection 

2. Hemisphere viewed from above the equator 
3. North polar projection 

4, South polar projection 


5. Hemisphere but viewed from some point other than the equator. 
Viewing position specified by XEAST and XWEST. 


An indicator which allows the user to reuse the mapping region 
defined in the previous call to the subroutine. IPLOT = 0 indicates 
that a new mapping region and/or a new scale is being defined with 
this call; therefore a new map is needed. IPLOT = 1 indicates that 
the same region as the previous call is being used, but that a new 
map of the region is desired. (This value of IPLOT should be used 
only after a call defining IEND = 0. See below.) Any other defini- 
tion of IPLOT indicates that this call defines new data which is to 
be plotted over the data plotted in the previous call. 


The plot advance control. IEND = O will plot the data for this call 
and advance the plotter to the origin of the next plot. Distance 
between plots = 4.0 inches. Any other definition of IEND will not 
terminate the plot and subsequent calls to the subroutine will cause 


the new data to be plotted over the older data. 


Contains the longitude axis label in Hollerith format. Unused if no 
label desired. 


Number of characters in IXAX. Must = 0 is no longitude label is 


desired or for polar projection. 


Contains the latitude axis label in Hollerith format. Unused if no 
label desired. 
Number of characters in IYAX. Must = 0 if no latitude label is 


desired for polar projection. 


Contains the plot title in Hollerith format. Unused if no label 
desired. 
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IT Number of characters in ITLE. Must = 0 if no plot title is desired. 
LABEL Contains any identification which the user wishes to apply to this 
set of data. The identifier must be in Hollerith format. It will 5 


be positioned above and along the line beginning with the first point 
inside the mapping region. - 


IL The number of characters in LABEL. Must = 0 if no data identifica- 
tion is desired. 


NOTE 1: The user must make provision for a physical TAPE3, EF20, 1/2 inch, 
556 BPI. 
NOTE 2: For the special case of using WORLD with Calcomp W2 paper, use the 


following parameters: XSIZE = 9.0, YSIZE = 18.0, XEAST = 360.0, 
XWEST = 0, YNORTH = 90.0, MAP = 1, IPLOT = 2. Other parameters may 
be set as appropriate. 
Considerable effort has been expended to make subroutine WORLD as general as 
the average programmer may need. If problems do occur or if help is needed on 
very specific problems, contact Captain Wiley, WLTH, AFWL. 
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SECTION IIT 


GENERAL COMMENTS TO THE USER 


The WORLD package consists of a main routine, WORLD, and 11 supporting 
routines. All control of the supporting routines is determined by the calling 
parameters of WORLD. WORLD has the following capabilities: 


1. To draw a map of any specified region and label it as desired. 

2. To draw a map and overlay data. 

3. To overlay several sets of data on same map. 

4. To ignore all data falling outside the specified map region. 

5. To do any of above in mercator, polar, or hemispherical projection. 


There are some limitations on the use of this routine. Although the calling 
sequence specifies both XSIZE and YSIZE, one or the other is always ignored. If 
YSIZE is 10.0 inches or less, then XSIZE will be adjusted to ensure that the 
latitude and longitude scaling factors are equal. If YSIZE is greater than 10.0 
inches, but XSIZE is 10.9 in-hes or less, then YSIZE will be adjusted to ensure 
that the scaling factors remain equal. (This will result in a reorientation of 
the axes and the values of XSIZE and YSIZE will be switched, which may cause 
problems in the calling program.) If both XSIZE and YSIZE are greater than 10.0 
inches, the size is incompatible with the Calcomp plotter; therefore WORLD will 
print an error message and stop execution. Even though the range of all longi- 
tudes can be -180 to +180 when calling WORLD, the returned values will be between 
0 and 360, as will be all longitude labels on the plot. In addition, the YLAT 
array will be altered when doing polar projections. When MAP = 5, XWEST, XEAST, 
XLONG, and YLAT are all changed to the new, rotated coordinate system, 


Features of mercator maps include a boundary around the map with tick marks 
every inch on all four sides, latitude values defined at each tick mark along the 
left edge with longitude values defined along the bottom, and any labels which 
the user provides. Polar projections will have the diameter of YSIZE. Quadrant 
longitudes will be drawn and labeled. North polar maps have 0° longitude at the 
bottom while south polar maps have 0° longitude at the top. Latitude circles 


will be drawn and labeled every inch from pole to boundary. 
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It must be pointed out to the user that each call to the routine which 
specifies that a new region is to be mapped involves rewinding a tape, reading 
that tape, and sorting the data. Therefore. one should use care in defining the 
calling perameters so that unnecessary tape handling can be avoided. Remember 
that one can change the scaling (i.e., the size) without doing any tape handling 
as long es the map region remains constant. (HINT: Copying the physical rape 
to a disk TAPE} will reduce the tape handling problen.) 


Bach of the major subroutines is briefly, described helow: 


SCALE Uses the calling parameters XEAST, XWEST, YNOKTH, YSOUTH end MAP 
to fie the map into the area specified by YSIZE and/or XSIZE and 
eet scaling factors for other subroutines. 


Plots the continenta) outliner ‘9 the :lssired projection. 


El!minates all continental areas outside of map region. 


Draws boundaries and labels plots. 
Treasforas the continental outline end plot data to the proper 
viewing engle. 


The normal starting point for all plots except Calcomp W2 is 1/2 inch fros 
the right edge of the Calcomp paper. For W2, the zero point is -90° latitude 
and 0° longitude (1.e., the crosshairs aust be at -90° latitude and 20° longitude). 


BRAY 

ELIM 

WLOT2 Scales and plots the user's dara over the world sap. 
TITLE 

VIEW 


Subroutine WORLD is used in conjunction with any of the Calcomp plot packages 
or eimulation routines. 


ERROR WECSAGTS 

SIZE 1S INCOMPATIBLE WITH CALCOMP PAPER SIZE. 

RBIZE © XX.X YSIZE © XX.xx 

eceure whenever both XSIZE and YSIZE are greater than 19.1 inches. 
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SAMPLE PLOTS 
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TEST OF MERCATOR 


bY 6s RD =F | 28 1 .e 188.6 


Figure 1. Modified Mercator 


CALL WORLD (DUM, DUM, 0., 6., 155., 355., 90., -90., 1, 0, 0, 
QULONGITUDE, 9, GULATITUDE, 8, 16NTEST OF MERCATOR, 16, DUM, 0) 
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TEST OF HEMISPHERE 


> a 


Figure 2. Hemisphere 


CALL WORLD (DUM, DUM, 0, 6., 6., 155., 355., 90., ~90., 2, 0, 0, 
DUM, 0, DUM, 0, l16HTEST OF HEMISPHERE, 18, DUM, 0) 
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WORLD (DUM, 


DUM, 0, 6., 6 


TEST OF NORTH, POLAR 


Figure 3. North Polar 


DUM, 0., DUM, 3, 0, 0, DUM, 


CALL » 6., 6., DUM, 
0, DUM, 0, 19HTEST OF NORTH POLAR, 19, DUM, 0) 
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TEST OF SGUTH POLAR 


=o 
180.0 


Figure 4. South Polar 


CALL WORLD (DUM, DUM, 0, 6., 6., DUM, DUM, DUM, 0., 4, 0, 0, DUM, 
0, DUM, 0, 19HTEST OF SOUTH POLAR, 19, DUM, 0) 


1l 
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TEST GF GFF-EQUATGR VIEW 


Figure 5. Hemisphere (Off-Equator View) 


CALL WORLD (DUM, DUM, 0, 6., 6., 250., 45., 90., -90., 5, 0, 0, 
DUM, 0, DUM, 0, 24HTEST OF OFF-EQUATOR VIEW, 24, DUM, 0) 
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- SATELLITE GROUND TRACK 


180.5 220.0 60.0 en.0 30.0 22.0 980.0 


Figure 6. Data Plot 
CALL WORLD (XLONG, YLAT, N, 6., 6., 360., 180, 90., -90., 1, 0, 0, 


QHLONGITUDE, 9, SHLATITUDE, 8, 22HSATELLITE GROUND TRACK, 22, 
6HPASS 1, 6) 
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